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fü 要 : 为 了 解 地 形 因子 对 树种 分 布 变化 的 影响 ,更 科学 地 进行 防护 林 空 间 配置 ,本 研究 利用 陕西 
省 延安 市 吴起 县 数字 高 程 模型 (DEM) 数 据 以 及 两 期 树种 分 布 数据 ,采用 Pearson 相关 分 析 和 随机 
森林 模型 ,探讨 陕西 省 延安 市 吴起 县 典型 造林 树种 [刺槐 (Robinia pseudoacacia ) | 1i Æ (Armeniaca si- 
birica)、 油 松 (Pinus tabuliformis ) fu W 3 (Hippophae rhamnoides)j] 与 地 形 因 子 ( 海 拔 、 坡 度 、 坡 向 、 曲 
率 、 地 形 起 伏 度 和 地 形 粗 糙 度 ) 的 关系 。 结 果 表 明 :(1) 2007 一 2019 年 刺槐 、 山 查 和 油 松 的 分 布 呈 增 
加 趋势 ,研究 区 内 分 别 在 83.373%、73.271% 和 72.415% 的 区 域 呈 增 加 趋势 ; 沙 琼 的 分 布 呈 减 少 趋 
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35 2€ 61.737908] X 39, 3E 38 , (2) 38 E Pearson 相关 分 析 和 随机 森林 重要 性 检验 ,海拔 和 坡 向 是 影响 
黄土 高 原 地 区 典型 造林 树种 分 布 的 主要 地 形 因 子 。(3) 4 种 主要 造林 乔 灌木 树种 刺槐 、 山 查 、 油 松 和 
沙 琼 适 宜生 长 海拔 范围 分 别 为 1410.105~1510.215 m.1444.873-1553.210 m,1380.132-1451.762 m 
和 1365.621-1562.348 mo (4) 刺槐 适宜 生长 坡 向 范围 为 0~20.520° 、50.280°~62.310° 、82.440?°~ 
95.280? fu 277.280? 360.000? ; 山 查 适宜 生长 坡 向 范围 为 0~81.270°? 和 265.310°~360.000?; 油 松 适宜 


生长 坡 向 范围 为 0~100.180°; 沙 球 适 宜生 长 坡 向 范围 为 0~74.540? 和 280.820°~360.000?°。 


关 键 词 : 树种 分 布 ; 地 形 因子 ; 随机 森林 模型 ; 黄土 高 原 


文章 编号 : 


森林 是 生态 系统 的 重要 组 成 部 分 ,在 拦截 降 
雨 .维持 径流 .调节 气候 防治 沙漠 化 和 水 土 保持 等 
方面 发 挥 重 要 作用 "…。 营 造 防 护林 逐渐 成 为 各 国 解 
决 生态 环境 问题 的 重要 手段 ,我 国 自 20 世 纪 50 代 
初 开 始 实 施 各 项 防护 林 体系 建设 工程 中 ,但 造林 初 
期 局 部 地 区 未 能 较 好 地 遵循 “ 适 地 适 树 ”原则 ,导致 
部 分 区 域 形成 的 防护 林 空间 配置 不 合理 。 因 此 ,为 
实现 防护 林 体系 建设 工程 的 可 持续 性 发 展 ,促进 防 
护林 结构 的 优化 ,进而 提升 防护 林 的 质量 ,探究 影 
响 防 护林 体系 营造 树种 分 布 的 因素 是 十 分 迫切 且 
有 意义 的 研究 工作 。 

目前 ,已 有 大 量 针 对 树种 分 布 与 其 影响 因子 之 
间 耦 合 关系 的 研究 成 果 ,认为 气候 土壤、 地 形 和 地 
质 等 “为 限制 树种 分 布 的 主要 因素 ,而 在 不 同 区 域 
尺度 范围 下 ,各 限制 因子 的 所 占 权 重 有 所 不 同 , 如 
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在 区 域 及 全 球 的 范围 内 ,气象 因子 是 影响 树种 分 布 
最 主要 的 因子 ,其 中 降水 .温度 是 主要 影响 因子 ”; 
地 形 因子 中 海拔 坡度 与 坡 向 是 主要 影响 因子 “”。 
现 有 研究 大 多 是 选取 样 地 点 进行 数据 调查 ,同时 考 
虑 气候、 地 形 和 土壤 等 因子 。 长 期 以 来 ,人 们 对 树 
种 分 布 驱 动因 素 的 研究 主要 侧重 于 气候 .土壤 因子 
与 地 形 因子 综合 作用 的 研究 ,而 局 域 范围 内 地 形 因 
子 更 能 反映 树种 的 生长 状况 ,选取 地 形 因子 为 指标 
指导 造林 计划 的 制订 是 最 快捷 ,简便 的 方法 。 现 有 
研究 对 象 大 多 选用 东北 .西南 和 青藏 高 原 等 地 区 的 
天 然 次 生 林 树 种 , 较 少 对 黄土 高 原 地 区 造林 树种 分 
布 的 驱动 因素 进行 研究 。 因 此 人 研究 黄土 高 原 地 区 
人 工 造林 局 域 范围 内 的 地 形 因 子 是 如 何 影响 树种 
分 布 的 问题 具有 重要 意义 。 

因此 ,为 探究 各 树种 在 不 同 地 形 因 子 下 的 分 布 
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规律 ,准确 理解 局 部 范围 内 地 形变 化 对 黄土 高 原 树 
种 分 布 的 影响 。 基 于 2007 年 和 2019 年 树种 分 布 数 
据 , 对 陕西 省 延安 市 吴起 县 防护 林 主 要 造林 树种 
[刺槐 (Robinia pseudoacacia ) , 1156 (Armeniaca sibiri- 
ca) PA (Pinus tabuliformis ) Vb (Hippophae rham- 
noide) | 的 分 布 进 行 了 时 空 变化 分 析 , 并 应 用 随机 和 森 
林 模 型 分 别 对 4 种 造林 树种 与 地 形 因 子 之 间 的 关系 
进行 定量 分 析 , 以 期 为 优化 空间 配置 ,合理 进行 造 
林 项 目 提供 理论 依据 。 


1 数据 与 方法 


1.1 研究 区 概况 

人 研究 区 位 于 水 蚀 风蚀 交错 区 的 陕西 省 延安 市 
吴起 县 (107"38'57"~108。32' 49" E, 36? 33' 33" ~ 
37?24'27"N) ,面积 98.980 km’, YFIR 1233-1809 mo 
研究 区 处 于 半 和 干旱 温带 大 陆 性 季风 气候 区 ,多 年 平 
均 降 水 量 为 466.700 mm ,平均 气温 为 8.050 *C ,年 均 
地 面 蒸发 量 为 891.200 mm。 自 1998 年 起 开始 退耕 
还 林 还 草 工程 以 来 ,研究 区 在 荒山 及 坡地 营造 水 土 
保持 林 , 主要 造林 树种 有 油 松 、 刺 槐 、 山 柑 URRE, 
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图 1 HUE 


典型 配置 方式 有 油 松 纯 林 、 刺 槐 纯 林 、 山 禁 纯 林 、 沙 
环 纯 林 V UP VP ROIG AO 、 刺 槐 沙棘 混交 林 和 山林 
沙棘 混交 林 。 为 科学 地 考虑 地 形 因 子 对 树种 分 布 
的 影响 ,选取 2002 年 以 前 完成 造林 且 受 人 为 因素 干 
扰 较 小 的 梁 坡 、 沟 坡 和 梁 踢 顶 自然 生长 的 防护 林 进 
行 研究 。 
1.2 数据 来 源 

利用 数字 高 程 模型 (DEM ) ,提取 研究 区 的 地 形 
因子 ” ,包括 :海拔 坡度 . 坡 向 、 曲 率 、 地 形 起 伏 度 
以 及 地 形 粗糙 度 (Topographic roughness index, TRI) 
(图 1)。 研 究 区 的 树种 分 布 信息 来 源 于 森林 资源 二 
类 调查 及 野外 调查 数据 (野外 调查 进行 时 间 为 2019 
年 7 一 9 月 )。 
1.3 树种 分 布 变化 等 级 划分 

为 对 树种 分 布 变 化 情况 进行 可 视 分 析 , 更 准确 
地 描述 树种 分 布 情况 ,结合 研究 区 内 实际 情况 , 根 
据 两 期 树种 分 布 信息 ,对 每 个 地 块 内 的 树种 进行 赋 
值 处 理 。 地 块 内 某 树 种 为 第 一 优势 树种 赋值 3 ,为 
第 二 优势 树种 赋值 2, 为 第 三 优势 树种 赋值 1。 分 别 
剔除 无 该 树种 分 布 区 域 ,以 2019 年 与 2007 年 的 差 
值 来 表示 树种 分 布 变化 情况 的 。 将 研究 区 树种 变 


因子 提取 结果 


Fig. 1 Extraction results of terrain factor 
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化 情况 划分 为 7 个 等 级 ,依次 为 :显著 减少 (-3) \ 减 
少 (-2) \ 稍 减少 (-1) 基本 不 变 (0) 、 稍 增加 (1) 385 
加 (2) 和 显著 增加 (3)。 
1.4 研究 方法 

本 研究 采用 随机 森林 模型 分 析 地 形 因子 与 树 
种 分 布 的 关系 。 随 机 森林 模型 是 由 Breiman" F 
2001 年 提出 的 一 种 统计 学 习 理论 ,该 方法 的 核心 思 
想 是 组 合 多 个 决策 树 ,对 变量 进行 判别 ,在 分 析 自 
变量 与 因 变 量 相关 性 时 精度 高 ,可 直接 估计 变量 重 
要 性 ,并 旦 在 建 模 过 程 中 产生 无 偏 估计 ,可 避免 过 
度 拟 合 。 同 时 ,该 模型 在 处 理 大 量 信息 时 ,即使 有 
部 分 数据 缺失 ,也 不 影响 其 建 模 精 度 ,在 处 理 大 数据 
集 上 优越 于 其 他 算法 ™, 是 一 种 常用 的 处 理 大 量 净 
据 和 多 个 变量 的 数学 方法 ,被 广泛 应 用 于 生态 学 ” 
等 领域 。 

随机 森林 模型 以 总 体 分 类 精度 和 Kappa 系数 作 
为 判别 模型 精确 度 和 可 靠 性 的 指标 。Kappa 系数 的 
判别 标准 如 表 1 所 示 。 本 人 研究 选取 随机 森林 模型 
对 各 树种 地 形 因 子 与 树种 分 布 进行 分 析 , 模 型 精度 
均 较 高 。4 种 树种 的 随机 森林 模型 总 体 精度 均 在 
75.000% 以 上 ,Kappa 系 数值 在 0.600 以 上 ,表明 选取 
的 地 形 因子 能 够 解释 影响 树种 分 布 变化 的 因素 ,应 
用 随机 森林 模型 分 析 树 种 分 布 与 地 形 因 子 的 关系 
结果 可 靠 。 


表 1 Kappa 系 数 与 精度 对 应 关系 
Tab.1 Correspondence between Kappa 


coefficient and accuracy 


Z2 不 同 mtry 参数 对 应 误差 


Tab.2 Corresponding errors of different mtry 


parameters 

特征 数 me — 
刺槐 iif 油 松 沙棘 

1 0.231 0.209 0.196 0.229 

2 0.218 0.200 0.177 0.219 

3 0.207 0.195 0.172 0.215 

4 0.208 0.193 0.169 0.214 

5 0.204 0.188 0.166 0.214 

6 0.203 0.186 0.170 0.213 


选 值 需 结合 模型 错误 率 与 决策 树 数 量 的 关系 。 
选取 错误 率 最 低 的 mtry 建 模 ,得 到 决策 树 数量 与 错 
误 率 关系 。 根 据 结果 最 终 确 定 刺 槐 以 mtry=6， 
ntree=3000; 山 #7 以 mtry=6, ntree=4100; 油 松 以 
mtry=5 ,ntree=2700; P ili A mtry=6 , ntree-2600 ,建立 
随机 和 森林 模型 。 

随机 森林 算法 中 的 相对 重要 性 (Mean decrease 
accuracy) 可 比较 变量 的 重要 性 ,节点 纯度 (Mean de- 
crease gini) 可 描述 变量 对 分 类 树 每 个 节点 上 观测 值 
异 质 性 的 影响 。 相 对 重要 性 和 节点 纯度 值 越 大 表 
示 该 变量 的 重要 性 越 大 。 因 此 ,本 研究 选取 随机 森 
林 算 法 中 的 相对 重要 性 指数 .节点 纯度 指数 来 确定 
各 地 形 因子 对 树种 分 布 变 化 的 影响 程度 。 


2 结果 与 分 析 


21 各 树种 分 布 变化 情况 分 析 
与 2007 年 相 比 ,2019 年 油 松 在 83.372% 的 区 域 


Kappa 系数 精度 评价 
0.800~1.000 几 近 完美 
0.600-0.800 精确 
0.400~0.600 较 好 
0.200~0.400 正常 
0.000~0.200 较 差 
<0.000 非常 差 


随机 森林 中 有 2 个 重要 参数 mtry 与 ntree。 本 人 研 
究 选取 6 个 自 变 量 , 因 此 以 1 为 单位 步 长 , 设 定 mtry 
参数 自 1-6 进行 逐步 长 遍历 ,选择 错误 率 最 低 的 
miy 值 为 参数 。 当 刺槐 A .以 及 沙棘 的 特征 数 为 
6 时 ,随机 森林 模型 的 错误 率 最 低 , 油 松 的 特征 数 为 
5 时 ,随机 森林 模型 的 错误 率 最 低 ( 表 2)。 

参数 ntree 是 建 模 时 决策 树 的 数量 ,过 低 会 导致 
模型 错误 率 偏 高 ,过 高 会 无 故 耗 费 模型 运行 时 间 ， 


呈 增 加 趋势 ( 稍 增加 .增加 和 显著 增加 ) ;其 次 为 刺 
EMIA ,分 别 在 73.272% 和 72.615% 的 区 域 呈 增加 
趋势 ;增加 最 少 的 为 沙 环 , 仅 在 23.555% 的 区 域 呈 增 
加 趋势 ( 表 3)。 沙 环 维 持 基本 不 变 的 比例 最 大 ,为 
14.709% ; IIT YE 12.481% 区 域 基本 不 变 ; 而 刺槐 和 
油 松 相 对 来 说 基本 不 变 区 域 较 少 ,分 别 为 9.887% 和 
9.59690, WIRA 61.738% H DX Tk, E Va / b ES A CRI UR 
少 减少 和 显著 减少 ) ,在 4 类 树种 中 衰退 趋势 最 为 
明显 ;其 次 为 刺槐 与 山 柑 ,分 别 在 16.842% 和 
15.104906 的 区 域 呈 减少 趋势 ;最 后 是 油 松 , 仅 在 
7.031% 的 区 域 呈 减少 趋势 。 

图 2 显示 了 研究 区 树种 分 布 空间 变化 情况 。 刺 
槐 在 整个 区 域内 分 布 范 围 较 广 , 且 几 种 变化 类 型 分 
布 区 域 较为 均匀 ,比例 相当 。 山 杏 广泛 且 均 匀 分 布 
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表 3 树种 分 布 面积 的 变化 


Tab.3 Changes in the distribution area of tree species 


m m" — 刺槐 — IA — 油 松 Vj 
面积 km ”百分比 /% 面积 km” ”百分比 /% j 积 /km ”百分比 /% 面积 km” ”百分比 /% 

-3 I 显著 减少 1.627 8.557 0.706 4.487 0.450 3.973 9.076 35.497 
-2 T ma 0.358 1.882 0.811 5.159 0.137 1.212 2.451 9.585 
-1 M 稍 减少 1.217 6.404 0.858 5.458 0.209 1.846 4.259 16.656 
0 IV 基本 不 变 1.879 9.887 1.963 12.481 1.086 9.596 3.761 14.709 
1 V 稍 增加 2.304 12.123 2.176 13.837 1.161 10.258 1.826 7.141 
2 VE 增加 3.941 20.730 4.636 29.482 2.224 19.648 2.034 7.957 
3 WL 显著 增加 7.683 40.419 4.576 29.096 6.051 53.466 2.162 8.457 


mm | uv» ru V 稍 增加 
mm D 减少 mm V 增加 uer 
mm QD Rb —— mE VI 显 著 增加 i 


[3 IV 基本 不 变 


图 2 树种 分 布 变化 


Fig.2 Changes oftree species distribution 


于 整个 研究 区 。 油 松 集中 分 布 于 人 研究 区 北部 ,其 他 
区 域 少 有 分 布 ,变化 类 型 丰 分 布 最 为 广泛 ,表明 油 
松 的 增加 趋势 最 明显 。4 类 树种 中 沙 环 分 布 最 广 ， 
几 种 变化 类 型 中 工 类 分 布 最 多 , 且 均 匀 分 布 于 整个 
研究 区 ,表明 沙 环 整体 衰退 。 


22 树种 分 布 变化 与 地 形 因子 的 相关 性 

对 研究 区 内 的 刺槐 山林 . 油 松 和 沙 环 与 6 种 地 
形 因 子 进 行 Pearson 相关 性 分 析 ,分 析 结 果 如 表 4 所 
示 。 和 刺槐 的 分 布 与 海拔 和 曲率 显著 正 相 关 ,相关 性 
系数 分 别 为 0.087 和 0.033。 山 查 的 分 布 与 海拔 显著 
正 相 关 , 相 关 性 系数 为 0.153, 与 起 伏 度 和 坡度 显著 
负 相 关 , 相 关 性 系数 分 别 为 -0.039 和 -0.035。 油 松 
的 分 布 与 坡度 .粗糙 度 和 起 伏 度 显著 正 相 关 ,相关 
系数 分 别 为 0.053 .0.052 和 0.048。 沙 环 的 分 布 与 海 
拔 和 曲率 显著 负 相 关 ,相关 系数 分 别 为 -0.031 和 
-0.048 ,与 坡 向 显著 正 相 关 ,相关 系数 为 0.047。 结 
RRI ,海拔 是 影响 刺槐 、 山 柑 和 沙 环 分布 最 重要 
的 地 形 因 子 。 


表 4 树种 分 布 变化 与 地 形 因子 相关 性 分 析 
Tab.4 Correlation analysis between tree species 


distribution and topographical factors 


树种 MAE ”海拔 曲率 坡度 ERRE Wh 
刺槐 -0.020 0.087” 0.033” -0.026 -0.028 -0.008 
Hifi -0.028 0.153" 0.011  -0.035' -0.039' 0.108" 


油 松 0.052" 0.002 
ybi 0.015  -0.03L 


注 :*、** 分 别 表示 在 0.05 .0.01 水 平 ( 双 侧 ) 上 显著 相关 。 


-0.039 0.053” 0.048 -0.003 
-0.048" 0.020 0.031 — 0.047" 


2.3 地 形 因 子 对 树种 分 布 的 影响 程度 分 析 

根据 有 无 各 树种 分 布 样 点 对 变量 的 重要 性 ( 表 
5) 可 知 ,在 刺槐 分 布 区 域 ,对 有 无 树木 分 布点 最 为 
重要 的 地 形 因子 均 为 海拔 。 在 山 查 和 油 松 分 布 区 
域 ,对 无 树木 分 布点 区 域 最 重要 的 地 形 因 子 为 海 
拨 , 对 有 树木 分 布点 区 域 最 重要 的 地 形 因子 均 为 坡 
向 。 在 沙 环 分 布 区 域 ,对 有 无 树木 分 布点 最 为 重要 
的 地 形 因子 为 坡 向 。 表 明海 拔 和 坡 向 是 影响 4 类 树 
种 分 布 最 为 重要 的 地 形 因 子 。 
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表 5 变量 对 不 同 区 域 样 点 的 重要 性 


Tab.5 Importance of variables to sample points in different regions 
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地 形 因子 刺槐 IA 油 松 Vj 
无 树木 分 布 ”有 树木 分 布 无 树木 分 布 ”有 树木 分 布 无 树木 分 布 ”有 树木 分 布 无 树木 分 布 ”有 树木 分 布 
海拔 148.400 197.420 150.380 177.820 133.810 64.572 33.054 175.550 
坡度 30.620 52.054 24.813 40.814 34.691 37.926 38.650 44.833 
坡 向 113.890 72.196 117.700 217.450 86.055 122.800 198.240 206.110 
曲率 5.973 44.247 4.350 34.044 13.792 27.561 6.939 54.817 
起 伏 度 11.449 54.782 19.990 42.130 34.848 32.115 40.107 47.462 
粗糙 度 38.066 67.327 35.679 42.943 33.769 35.709 54.792 58.836 
注 :将 2007 年 为 正 值 但 2019 年 值 为 0 的 点 定义 为 无 某 树种 分 布 ,将 2007 年 和 2019 年 均 为 正 值 的 点 定义 为 有 某 树种 分 布 。 


由 各 树种 地 形 因 子 的 相对 重要 性 和 节点 纯度 
指数 排序 (图 3) 可 见 , 刺 槐 、 山 杏 ` 油 松 和 沙棘 2 个 重 


的 2 个 因子 绘制 地 形 因 子 对 各 树种 分 布 的 影响 图 
(图 4、 图 5), 纵 坐标 表示 地 形 因 子 对 树种 分 布 的 影 


要 性 指数 较 大 的 地 形 因子 均 是 海拔 和 坡 向 。 除 山 
查 的 节点 纯度 重要 性 指数 排名 最 后 的 为 坡度 ,其 余 
树种 重要 性 指数 排名 最 后 的 均 为 曲率 。 

计算 变量 对 有 无 树种 分 布 样 点 的 重要 性 、 相 对 


响 程 度 , 值 大 于 0 表示 为 适宜 分 布 , 值 小 于 0 表示 不 
适宜 分 布 ,影响 程度 随 绝对 值 的 增加 而 增加 。 

各 树种 海拔 与 其 影响 程度 的 关系 整体 均 呈 现 
先 增加 后 减少 的 趋势 ,表明 各 树种 的 海拔 对 树种 分 


重要 性 和 节点 纯度 的 平均 值 , 可 得 到 地 形 因子 的 重 
要 性 排序 ,对 刺槐 影响 由 大 到 小 分 别 为 海拔 > 坡 向 > 
粗粮 度 > 起 伏 度 > 坡度 > 曲率 ;对 山 杏 影响 由 大 到 小 
分 别 为 坡 向 > 海拔 > 粗糙 度 > 起 伏 度 > 坡度 > 曲率 ;对 


布 的 影响 均 随 着 海拔 的 增加 而 增 大 ,到达 极 值 后 随 
着 海拔 的 增加 而 减 小 (图 4)。 刺 槐 在 海拔 1410.105 
m 以 下 和 1510.215 m 以 上 的 受 影响 程度 均 为 负 , 在 
海拔 1410.105~1510.215 m 受 影响 程度 为 正 ;山林 在 


油 松 影响 由 大 到 小 分 别 为 坡 向 > 海拔 > 坡度 > 起 伏 度 
> 粗糙 度 > 曲率 ;对 沙棘 影响 由 大 到 小 分 别 为 海拔 > 
坡 向 > 起 伏 度 > 粗糙 度 > 坡度 > 曲率 。 选 取 排名 靠 前 


(a1) 刺槐 相对 重要 性 (a2) 刺槐 节点 纯度 
曲率 。 曲率 [， 
pis 坡度 | ， 
th ip 
1l 坡度 ES 粗糙 度 [ 
坡 向 坡 向 上 
海拔 海拔 上 
0 100 200 300 100 200 300 400 500 
相对 重要 性 节点 纯度 
(c1) 油 松 相对 重要 性 (c2) 油 松 节 点 纯度 
曲率 曲率 。 
坡度 上 
人 次 
Tasa] 
xi 粗糙 度 B 
海拔 上 
坡 向 上 


0 50 100 150 200 


相对 重要 性 节点 纯度 


0 50 100 150 200 


海拔 1445.873 m 以 下 和 1553.210 m 以 上 受 影 响 程度 
均 为 负 ,在 海拔 1444.873~1553.210 m 受 影响 程度 为 
正 ; 油 松 在 海拔 1380.132 m 以 下 和 1451.762 m 以 上 


(b1) 山 杏 相 对 重要 性 (b2) 山 杏 节 点 纯度 
[UE 坡度 i 
坡度 曲率 
Mb sun bs auk 
R 起 伏 度 g 起 伏 度 
四 粗粮 度 3 粗糙 度 
海拔 海拔 
坡 向 坡 向 
0 100 200 300 0 100 200 300 
相对 重要 性 节点 纯度 
(d1) 沙棘 相对 重要 性 (d2) 沙 环 节点 纯度 
曲率 [「 曲率 『 
坡度 上 ， 
Ms ad 
A 起 伏 度 
粗糙 度 
坡 向 
海拔 


0 100 200 300 
相对 重要 性 


100 200 300 400 500 
节点 纯度 


图 3 地 形 因子 重要 性 指数 排序 
Fig. 3 Importance index ranking of topographical factors 
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图 4 海拔 对 各 树种 的 影响 
Fig.4 Effect ofaltitude on tree Species 
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图 5 坡 向 对 各 树种 的 影响 


Fig. 5 Effect of aspect on tree species 


受 影响 程度 为 负 ,在 海拔 1380.132~1451.762 m 受 影 
响 程 度 为 正 ; 沙 环 在 海拔 1365.621 m 以 下 和 
1562.348 m 以 上 受 影响 程度 为 负 ,在 海拔 1365.621~ 
1562.348 m 受 影响 程度 为 正 。 

各 树种 坡 向 与 其 影响 程度 的 关系 整体 均 呈 现 


Ji E HE EJ 74.5409 280.8209 22 E np FE BE H fa , 其 他 
坡 向 受 影响 程度 为 正 。 


3 讨论 


地 形 是 导致 各 种 生态 现象 和 过 程 发 生 分 异 的 
重要 因素 ,海拔 .坡度 \ 坡 向 等 各 因子 通过 对 地 表 物 
质 与 能 量 的 再 分 配 共同 决定 着 地 表 植 被 的 分 布 格 
局 “5 。 本 研究 利用 随机 森林 模型 评估 了 各 地 形 因 
子 对 树种 分 布 的 重要 性 ,结果 表明 海拔 和 坡 向 是 影 
响 刺 槐 , LLL ` 油 松 和 沙棘 分 布 的 主要 地 形 因子 。 
这 是 由 于 海拔 和 坡 向 导致 了 太阳 辐射 与 蒸 散 发 的 
不 同 ,进而 对 区 域 间 的 土壤 水 分 状况 造成 影响 ,而 
土壤 水 分 又 是 该 地 区 植物 生长 发 育 的 主要 限制 因 
子 呈 ,从 而 导致 了 生物 个 体 的 生长 发 育 与 树种 分 布 
的 不 同 ,使 得 海拔 和 坡 向 成 为 限制 研究 区 域 典 
型 造林 树种 分 布 的 关键 因子 。 

海拔 对 黄土 高 原 区 域 典 型 树种 分 布 的 影响 最 
大 ,这 与 谢 玉 彬 等 所 在 探究 选取 地 形 因子 中 的 海 
拔 .坡度 和 凸 度 对 天 童 地 区 常 绿 树种 和 落叶 树种 分 

影响 的 研究 结果 相 一 致 。 具 体 来 说 各 树种 适宜 生 
长 海拔 范围 由 大 到 小 依次 为 沙棘 (196.727 m) 山林 
(108.337 m) HJE (100.110 m) 和 油 松 (71.630 m) 。 
JK Ie TRES" BG AECOU SK ESECU BREUI 
研究 结果 表明 ,在 黄土 高 原 地 区 , 阴 坡 相 较 于 阳 坡 
更 适宜 植被 的 生长 ,本 研究 的 研究 结果 与 其 相 一 
致 , 即 4 类 树种 在 阴 坡 分 布 均 优 于 阳 坡 。 同 时 , 现 有 
人 研究 大 多 将 坡 向 分 为 阴 坡 (0~45°* 和 315°~360°)、 半 
阴 坡 (45°~135°) 、 阳 坡 (135°~225°) 、 半 阳 坡 (225°~ 
360° ) 等 级 ,选取 位 于 不 同 坡 向 等 级 样 地 作为 代表 进 
行 研 究 , 表 明 处 于 阴 坡 和 半 阴 坡 的 油 松 更 新 能 力 明 
显 优 于 阳 坡 ?3 , 沙 环 的 生长 状况 由 好 到 坏 顺 序 为 
阴 坡 、 半 阴 坡 、 半 阳 坡 、 阳 坡 ,这 与 本 人 研究 的 结果 
基本 一 致 , 即 造 林 时 树种 应 主要 种 植 于 水 分 条 件 较 
好 北 坡 , 均 不 适宜 种 植 在 南 坡 , 且 沙棘 和 山 奋 适宜 
生长 坡 向 范围 较 广 于 刺槐 和 油 松 。 此 外 ,本 研究 给 


先 减少 后 增加 的 趋势 ,表明 各 树种 的 坡 向 对 树种 分 
布 的 影响 均 随 着 坡 向 的 增加 而 减 小 ,到 达 极 值 后 随 
着 坡 向 的 增加 而 增加 (图 5)。 刺 槐 在 坡 向 20.520?~ 
50.280? 、62.310°~82.440° .95.280°~277.280° 受 影响 
程度 为 负 , 其 他 坡 向 为 正 ; 山 查 在 坡 向 81.270°~ 
265.310° 受 影响 程度 为 负 , 其 他 坡 向 为 正 ; 油 松 在 坡 
向 100.180? 以 上 受 影响 程度 为 正 , 其 他 坡 向 为 负 ; 沙 


出 刺槐 、 山 查 、 油 松 和 沙棘 适宜 生长 坡 向 范围 的 具 
体 数值 ,可 为 防护 林 工 程 后 续 建设 过 程 中 的 林 分 空 
间 结 构 优 化 和 可 持续 性 发 展 提 供 理论 依据 。 

对 研究 区 典型 造林 树种 的 时 空 分 布 进行 分 析 ， 
结果 表明 刺槐 . 山 杏 以 及 油 松 都 呈 增 加 趋势 ,只 有 
沙 束 呈 衰 退 趋 势 , 这 与 张 恰 只 等 人 研究 结果 一 致 。 
通过 随机 森林 模型 结果 显示 适宜 沙 环 生长 的 坡 向 
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为 阴 坡 ,海拔 为 1365.621~1562.348 m, SE Scb lE 
退 的 原因 可 能 是 在 半 干 旱 的 黄土 高 原 地 区 ,土壤 水 分 
限制 了 各 种 植被 的 生长 发 育 , 该 地 区 不 同 地 形 部 位 的 
土壤 水 分 含量 差异 较 大 ,在 初期 造林 时 ,作为 先锋 树 
种 的 沙 环 被 种 植 于 半 干 旱地 区 的 各 海拔 与 坡 向 ,在 不 
适宜 沙 为 生长 区 域 种 植 该 物种 导致 沙 环 过 度 利 用 水 
资源 而 导致 失去 再 生 能 力 甚 至 大 片 死亡 * > 。 


4 结论 


对 陕西 延安 市 吴起 县 约 100 km 范围 的 树种 分 
布 时 空 变化 特征 进行 了 研究 ,在 区 域 尺 度 上 对 地 形 
因子 与 树种 分 布 关系 进 行 定量 研究 ,获得 黄土 高 原 
常见 造林 树种 刺槐 LA ` 油 松 和 沙 环 适宜 生长 的 海 
拨 与 坡 向 具体 数值 范围 ,为 促进 森林 结构 的 优化 和 
提高 森林 质量 提供 理论 依据 。 主 要 结论 如 下 : 

(1) 与 2007 年 树种 分 布 相 比 , 油 松 分 布 面积 增 
加 最 多 , 旦 增加 区 域 主要 分 布 在 研究 区 北部 ;而 沙 
束 分 布 整体 呈 衰 退 趋势 ,退化 情况 呈现 于 整个 区 
域 , 曾 经 作为 先锋 树种 的 沙 环 逐渐 被 其 他 乔木 树种 
代替 ,表明 经 过 十 几 年 的 恢复 研究 区 植被 空间 结构 
逐步 合理 。 

(2) 对 刺槐 LA ` 油 松 和 沙 环 分 布 最 重要 的 地 
形 因子 均 为 海拔 和 坡 向 。 因 此 ,在 植树 造林 时 ,应 
根据 不 同 海拔 和 坡 向 选择 不 同 造林 树种 。 

(3) 刺槐 在 海拔 1410.105~1510.215 m 适宜 分 
布 , 山 柑 在 海拔 1444.873~1553.210 mn 适宜 分 布 , 油 
松 在 海拔 1380.132~1451.762 mi tr^) , Vb RETE TE 
4& 1365.621-1562.348 mm 适宜 分 布 。 各 树种 适宜 生 
长 海拔 范围 由 大 到 小 依次 为 沙棘 (196.727 m) MUA 
(108.337 m) 刺槐 (100.110 m) 和 油 松 (71.630 m) 。 

(4) 刺槐 适宜 分 布 坡 向 范围 为 0~20.520°、 
50.280°~62.310° 、82.440°~95.280° 和 277.280°~ 
360.000? ; 山 查 适宜 分 布 坡 向 范围 为 0~81.270° 和 
265.310°~360.000° ; 油 松 适宜 分 布 坡 回 范 围 为 0~ 
100.180? ; 1» R x& "EC 4) p 3E e 6 Fe] y 0~74.540° ffl 
280.820?~360.000"。 各 树种 适宜 生长 范围 由 大 到 小 
H ip d (180.440? ) , vb d (153.8209 ) 、 刺 槐 
(105.430?) Ih (100.180?) 。 
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Effects of topographical factors on tree species distribution of shelter forest in 
Loess hilly region of northern Shaanxi 
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Abstract: As an important part of various ecosystems, forests play an important role in intercepting rainfall, 
maintaining runoff, regulating climate, preventing desertification, and soil and water conservation. The 
construction of shelterbelts has also gradually become an important means of solving the ecological and 
environmental problems in several countries. However, because of the excessive pursuit of increased green area 
in the early stages of afforestation and the failure to follow the "suitable site and suitable tree" principle, the 
quality of protective forests formed in the early stages of afforestation in China has not been high overall. Thus, it 
is of great significance to study how topographic factors affect the distribution of tree species to realize 
sustainable development in shelterbelt system construction projects, promote the optimization of shelterbelt 
structures, and improve the quality of shelterbelts. In this study, digital elevation model data, forest resource 
survey results, and tree species distribution data from Wuqi County, Shaanxi Province, China were collected to 
investigate the influence of topographical factors on vegetation distribution and promote ecological restoration 
and shelter forest construction. The relationships between typical afforestation tree species (Pinus tabuliformis, 
Robinia pseudoacacia, Armeniaca sibiric, and Hippophae rhamnoides) and topographic factors (elevation, 
aspect, slope, curvature, recursive damped least squares, and topographic roughness index) were analyzed using 
Pearson analysis and a random forest model. The results demonstrated the following: (1) The distribution area of 
R. pseudoacacia, A. sibiric, and P. tabuliformis exhibited an upward trend from 2007 to 2019, increasing by 
83.373906, 73.271% , and 72.415% , respectively. The distribution area of H. rhamnoides exhibited a declining 
trend, decreasing by 61.737975. (2) On the local scale, the influence of terrain factors on the distribution of tree 
species cannot be ignored, with the most important being altitude and aspect. Hence, different afforestation tree 
species should be selected according to altitudes and directions. (3) Specifically, the altitude ranges suitable for R. 
pseudoacacia, A. sibiric, P. tabuliformis, and H. rhamnoides were found to be 1410.105-1510.215 m, 1444.873- 
1553.210 m, 1380.132-1451.762 m, and 1365.621-1562.348 m, respectively. The aspects suitable for R. 
pseudoacacia include 0-20.520?, 50.280?—62.310?, 82.440?—95.280?, and 277.280?—360.000?, whereas those for 
A. sibiric are 0—81.270? and 265.310?—-360.000?. The suitable range for P. tabulaeformis is 0—100.180?, whereas 
those for H. rhamnoides are 0—74.540? and 280.820?—360.000*. 


Key words: distribution of tree species; topographical factors; random forest model; Loess Plateau 


